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The goal of this problem was to provide a detailed study of the accuracy 
of boundary treatments with a range of incidence angles including shear and a 
sonic point. 

There are three parts. In each we solve the linearized Euler equations on a 
prescibed domain: (—2, 2) x (0, 1) with initial conditions consisting of a pressure 
dipole, entropy and vorticity disturbances. Here aq = ±2 are the artificial 
boundaries, the speed of sound is scaled to 1, and we solve up to t = 64. 

PART 1 

For part 1 the base flow is a uniform subsonic flow skew to the boundaries: 

Pi = 0.3, U 2 = 0.4. (1) 

In addition, periodic boundary conditions are prescribed a; 2 . 

The exact solution is given by the following formulas: 


p = P(x i — U\t,x 2 — U 2 t,t), p = D{ aq - U\t 1 X 2 — U 2 t,t), 
u = U (aq — U\t, x 2 — U 2 t, <), v = V ( aq — U\t, x 2 — U 2 t , t), 


where 


P(aq,x 2 ,f) = ^Bi ^2 
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V(xi,x 2 ,t) = - / ^ — (xi,x 2 ,s)ds + V 0 (x 1 ,x 2 ), 
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U 0 (xi,x 2 ) = /i(aq) J -(z,x 2 ,0)dz-J —(z,x 2 ,0)dz, 
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Vq(xi,x 2 ) = —f[(xi) J J — (z,w,0)dzdw. 
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and 

r ik = ( *1 ~ *1 ,i ) 2 + (*2 - x 2 ,ik) 2 , 
r Sk = (*i ~ *i,s) 2 + (*2 - x 2 ,sk) 2 , 

( 0, xi < —1.9 

/l(*l) = < 1 — e -((®i + l-9)/2.5) 8 ) | a ; 1 | < IJ 

\ 1, xi > 1.9. 

The parameters Bi, fii, #i,j, x 2 t ik, S , ns, x i,S , x 2 ,sk are chosen so that, to 
a high degree of accuracy (11 digits), the initial data is supported on (—2,2) 
and the boundary conditions are satisfied. The integrals are evaluated using 
a combination of Gaussian quadrature and endpoint corrected trapezoid rules, 
again to high accuracy. The infinite sums are truncated after the point where 
their contributions are below machine precision. We also note that the jump in 
/i is approximately 4 x 10 -13 . 

Precisely we chose a dipole-like initial configuration for the pressure pulse: 
n = r 2 = —.95, Hi = M2 = 30, B 2 = -Bi = 1, 

* 1,1 = — * 1,2 = 0 . 1 , * 2,10 = * 2,20 = 1 / 2 . 
and for the entropy pulse: 


Ms = 12, S= 1, *i„s = 0, *2, so = 1/2. 

To guarantee periodicity we have: 

*2 ,ik, * 2 ,Sk = ^ + -00 < k < 00. 

We note that similar solutions have been used to test boundary conditions 
for the linearized Euler equations in [2] and for the scalar wave equation in [1]. 

PARTS 2 AND 3 

In part 2 the base flow is given by the subsonic Couette flow: 

U\ = Mx 2 , M = 0.9, U 2 = 0, (2) 

and in part 3 by the transonic Couette flow: 

U l = Mx 2 , M = 1.2, U 2 = 0. (3) 

For these problems we replace the periodic boundary conditions by the wall 
boundary condition, v = 0. The initial conditions are defined by the same 
functions and parameters as part 1 except that the image source locations x 2 t ik 
are determined to guarantee compatibility with the wall conditions. For k > 0: 

*2,i,fc+l — 2 *2,?, — k, *2,i, — (fc+l) *2 ,ik, 


NAS A/CP— 2004-2 12954 


36 



X2,S,k+l = 2 — X2 ,S,-k, *2,S,-(fc+l) = — ^2 ,Sfe- 

In this case we don’t have a code which evaluates an exact solution. Instead 
we use a well-resolved numerical solution on a sufficiently long domain to elimi- 
nate the influence of the boundaries. Tha basic numerical scheme is identical to 
the one we used to solve these and other benchmark problems, and is described 
in more detail elsewhere in the proceedings. In time we use a standard 4th order 
Runge-Kutta method with time step dt = 1/2000: 128,000 steps for the entire 
solution. Space derivatives are calculated using an 8th order difference scheme 
on a square grid with an extra point near the boundaries (added for stability). 
Thus the mesh in the domain [-L, L] x [0, 1] has (nx + 3) * ( ny + 3) points 
(nx = 128 * 2 * L, ny = 128; h x = h y = 1/128) 

The length of the domain is chosen so that reflection from the left and right 
boundaries causing possible errors would not come before time t = 64.0 

— ^ + (L - 2) > 64.0 

Hence, L = 44 for Problem 2 (M = .9), and L = 47 for Problem 3 (M = 1.2). 
We note that this required 385, 120 points in the transonic case. We have not 
fully assessed the accuracy of this solution, but preliminary comparisons with 
coarser mesh solutions suggests that it is accurate to more than three digits. 
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